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Abstract. The modified 



V 2r / n + r 

introduced by D. Zagier in 1998 are extended to the polynomial case by replac- 
ing B r by the Bernoulli polynomials B r (x). Properties of these new polyno- 
mials are established using the umbral method as well as classical techniques. 
The values of x that yield periodic subsequences £?2 n+1 (a;) are classified. The 
strange 6-periodicity of Sln+l' es t aDusrle d by Zagier, is explained by exhibit- 
ing a decomposition of this sequence as the sum of two parts with periods 
2 and 3, respectively. Similar results for modifications of Euler numbers are 
stated. 



1. Introduction 
The Bernoulli numbers B n , defined by the generating function 

e 1 — 1 i n\ 

are rational numbers with B 2n +i = for n > 1 and Bi — — \- The sequence {B n } 
has remarkable properties and it appears in a variety of mathematical problems. 
Examples of such include the fact that the Riemann zeta function 



(1.2) coo 



' — ' n° 

n=l 



evaluated at an even positive integer s = In is a rational multiple of 7r 2n , the factor 
being 

(L3) -j*r - T^yr (- 1 ) 

Their denominators are completely determined by the von Staudt-Clausen theorem: 
the denominator of B 2n is the product of all primes p such that p — 1 divides 2n 
(see [11] for an elementary proof). It is often the numerators of B 2n that are 
the objects of interest. It is a remarkable mystery that there is no elementary 
formula associated to them. These numerators appear in connection to Fermat's 
last theorem (see [14] ) and also in relation to the group of smooth structures on 
n-spheres (see [9], page 530 and [10] for details). 
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D. Zagier |21] introduced the modified Bernoulli numbers 

(1-4) K = t( n t r )^ n>0 

and established the following amusing variant of i?2n+i = 0: 

Theorem 1.1. TTie sequence Bln+i * s periodic of period 6 u>i£/i values 

r3 _i _1 1 1 _3-i 

t 4 i 41 4'4'4' 4 J ' 

One of the goals of this work is to extend this result to the polynomial 

/ s v—^ (n + r\ B r (x) 
V 2r J n + r 

which is defined here as the Zagier polynomial. Here B r (x) is the classical Bernoulli 
polynomial defined by the generating function 

n=0 

The objective of the paper is to produce analogues of standard results for B n {x) 
for the Zagier polynomials _B*(x). For example, a generating function for these 
polynomials appears in Theorem 15.11 as 

00 

(1.7) B n(*)* n =-^ogz--^(z + l/z-l-x) 

n=l 

where VK 2 -) = ^'( x )/^( x ) is t ne digamma function. The generating function (II. 7[) 
really corresponds to the less elementary expression 

00 



(1.8) VB n (z)z n = ic(2,l/«-a! + l) 

n=0 

Here C( s j°) is t ne Hurwitz zeta function, defined by 



00 

(I- 9 ) CM = £7—- -T7- 

A derivation of (|1.8[) is given in Sectional The next example corresponds to the 
derivative rule B' n (x) — nB n -i(x) for the Bernoulli polynomials. It appears in 
Theorem |Q as 

LfJ 

— - J2(2j - VB^x) for n even 



and 



, LfJ 

—B* n {x) = - + Y J 2jB* 2j (x) for n odd. 



Finally, the analogue of the symmetry relation B n (l — x) = (~l) n B n (x) is estab- 
lished as 

(1.10) B* n (-x- 3) = (-!)* B*(x). 

This is the content of Theorem 111.11 
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The original motivating factor for this work was to extend Theorem 1 1.1 1 to other 
values of -B^n+if^)- The main result presented here is a classification of the values 
x G K for which B 2n . i(x) is a periodic sequence, Zagier's case being x = 0. 

Theorem 1.2. Assume {-B^n+il 2 -)} is aperiodic sequence. Thenx £ {—3, —2, —1, 0} 



In the case of even degree, the natural result is expressed in terms of the difference 



Theorem 1.3. Assume {AZ, n (x)} is a periodic sequence. Then x £ {—1, 0, 1, 2}. 
The period is 3 for x = — 1 and x = 2 while A^x) vanishes identically for x = 
and x = 1. 

An outline of the paper is given next. Section [2] contains a basic introduction 
to the umbral calculus with a special emphasis on the operational rules for the 
Bernoulli umbra. Section [3] gives the generating function of the modified Bernoulli 
numbers i3* and this is used to give a proof of the 6-periodicity of -Bln+i- An 
inversion formula expressing B n in terms of B* is given in Section 2J The proof 
extends to the polynomial case. The generating function for the Zagier polynomial 
B*(x) is established in Section [5] and an introduction to the arithmetic properties 
of special values of these polynomials appears in Section [5] Expressions for the 
derivatives of the Zagier polynomials are given in Section [5J Some binomials sums 
employed in the proof of these results are given in Section [7J The basic properties 
of Chebyshev polynomials are reviewed in Section [S] and used in Section ITTJI to give 
a representation of the Zagier polynomials in terms of Bernoulli and Chebyshev 
polynomials and also to prove a symmetry property of B*(x) in Section [TT] Section 
1121 contains one of the main results: the classification of all periodic sequences of 
the form f?2 n +i(a0- This result extends the original theorem of D. Zagier on the 
6-periodicity of -B^n+i- Several additional properties of the Zagier polynomials are 
stated in Section [T3] The results discussed in the present paper can be extended 
without difficulty to the case of Euler polynomials. These extensions are stated 
in Section [14] and used in Section [15] to establish a duplication formula for Zagier 
polynomials. 



In the classical umbral calculus, as introduced by J. Blissard [T], the terms a n of 
a sequence are formally replaced by the powers o™ of a new variable a, named the 
umbra of the sequence {a n }. The original sequence is recovered by the evaluation 
map 



The introduction of an umbra for {a n } requires a constitutive equation that reflects 
the properties of the original sequence. These ideas are illustrated with the umbra 
23 of the Bernoulli numbers {B n }. 

An alternative approach to (|1.1[) . as a definition for the Bernoulli numbers B n , 
is to use the equivalent recursion formula 




Al n {x) = B* 2n {x) - B* 2n {-1). 



2. The umbral calculus 



(2.1) 



eval{a"} = a, 



(2.2) 




for n > 1, 



4 



A. DIXIT, V. H. MOLL, AND C. VIGNAT 



complemented by the initial condition Bo = 1. In terms of the Bernoulli umbra 03, 
this recursion is written as 

(2.3) -05 = 05 + 1. 

This is a constitutive equation for the Bernoulli umbra and the numbers B n are 
then obtained via the evaluation map 

(2.4) eval{03' 1 } = B n , 

The umbral method is illustrated by computing the first few values of B n , starting 
with the initial condition B = 1. The choice n = 2 in (|2.3p gives 

(2.5) 03 2 = (03 + l) 2 = 03 2 + 205 1 + 03°. 

The evaluation map then gives B2 = B2 + 2B\ + Bq that simplifies to 2B\ + Bq = 
and this yields B\ = —1/2. Similarly, n = 3 gives 

(2.6) OS 3 = (03 + 1) 3 = 03 3 + 303 2 + 303 1 + 05°, 

and the evaluation map produces £>3+3£>2 + 3£?i+-E>o = -B3 an d B% = 1 is obtained. 
The reader will find more details about these ideas in [BJ. 

The evaluation map of the Bernoulli umbra 03 may be defined at the level of 
generating functions by 

(2.7) eval{exp(<03)} = - 



1 

Similarly, the extension of (|2.7[) to the umbrae 03 (x) for the Bernoulli polynomials 
in (|1.6p is defined by 

(2.8) eval |exp(i03(x))} = —7- — — . 

e l — 1 

It is a general statement about umbral calculus that the operation eval is linear. 
Moreover, expressions independent of the corresponding umbra are to be treated 
as constant with respect to eval. Some further operational rules, particular for the 
Bernoulli umbra, are stated next. 

Lemma 2.1. The relation 

(2.9) eval{03(a;)} = eval{x + 03}. 
holds. 

Proof. This comes directly from 

te xt t 

(2.10) eval{exp(£03(x))} = and eval {exp(<05)l = . 

e l — 1 e 1 — 1 

□ 

Lemma 2.2. Let P be a polynomial. Then 

(2.11) eval {P{x + 03 + 1)} = eval {P(x + 03)} + P'{x). 

Proof. This is verified first for monomials using (|2.3p and then extended by linearity 
to the polynomial case. □ 
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The next step is to present a procedure to evaluate nonlinear functions of the 
Bernoulli polynomials. This can be done using the umbral approach but we intro- 
duce here an equivalent probabilistic formalism that is easier to use in a variety of 
examples. These two approaches will be described further in [5] where some of the 
results presented in [20 will be established using this formalism. 

The notation 

(2.12) E[h(X)]= [ h[x)f x (x)dx 

Jw 

is used here for the expectation operator based on the random variable X with 
probability density fx- The class of admissible functions h is restricted by the 
existence of the integral (|2.12l) . The equation (|2.13p shows that a probabilistic 
equivalent of the eval operator of umbral calculus is the expectation operator with 
respect to the probability distribution (|2 . 14[) . 

Theorem 2.3. There exists a real valued random variable Lb with probability 
density fi, B {x) such that, for all admissible functions h, 

(2.13) evaX{h(<B{x))} = E [h(x - 1/2 + iL B )] 

where the expectation is defined in (|2.12[) . The density of Lb is given by 

(2.14) f hB {x) = |sech 2 (7ra;), for x e R. 
Proof. Put x — in the special case 

(2.15) eval{exp(£Q3(x))} = E [exp (t(x - \ + iL B ))] 
of (|2.13p and use (|2 .8[) to produce 

(2.16) E [exp(itL B )] = ^csch Q 

Let /l b (x) be the density of Lb and write (|2.16[) as 

f°° t ft 

(2.17) J cos(tu)f LB (u)du = -csch ( - 

assuming the symmetry of Lb- The result now follows from entry 3.982.1 in [7] 

nt , / 7rt " 



f nt 
(2.18) / sech 2 (au) cos(tu) du = ^-csch 

J-oc a 



2a 



□ 

Note 2.4. The integral representation of the Bernoulli polynomials 
(2.19) B n (x) =E[(x- \+iL B ) n ] - 

is a special case of Theorem 12.31 The formula (|2 . 19[) is stated in non-probabilistic 
language as 



(2.20) B n (x) = - I (x-\+it) n sech 2 (7ri) dt. 



To the best of our knowledge, this evaluation first appeared in [T5]. The role played 
by sech 2 x as a solitary wave for the Kortweg-de Vries equation has prompted the 
titles of the evaluations of ([2^20]) in [2j [8] - 
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The next result uses Theorem l2.3l to evaluate a nonlinear function of the Bernoulli 
polynomials that will be needed later. More examples will appear in the companion 
paper [5]. 



Theorem 2.5. Letip{x) — T'(x)/T(x) be the digamma function. Then 

(2.21) eval{log«8(a;)} = V>(i + |a;-!|). 
In particular, for x > h , 

(2.22) eval {log 93(2;)} = i/)(x). 
Proof. Theorem 12.31 with h[x) — log a; gives 

(2.23) eval {log <B(x)} = E [log (x - \ + iL B )\ ■ 

The density f^ B is an even function, therefore the random variables Lb and — Lb 
have the same distribution. This gives 

eval {log <B(x)} = [log ((a 

= log(x-±) + h 
A linear scaling of entry 4.373.4 in [7] gives 

C TT\fb 



1)2 

2> 



log I I 



T 2 

B 



1)2 



(2.24) 



log(l + 6u 2 ) J 2 

; j d,U = - 

sinh cu c 



log 



for b, c > 0. Define 

(2.25) h(b,c) 

and observe that 

r °° log(l + bu 2 ) 



log 



TT\fb 



TT\fb 

■K\fb 
' ~~2c 



2c 



7TV& 







sinh 2ttu 



It follows that 
(2.26) 



du = — 
4 

1 

4 



log(l + bu 2 



irVb 

log(l + bu 2 ) du 



cosh 2 iiu sinh 2 iru 
log(l + bu 2 ) ( cosh 2 7ru 



COSh ITU 



sinh ttu 



— 1 ) du 



■ du = h(b, 7r) - 4/i(6, 27r). 
cosh 7ru 



Now take 6 = (x — |) 2 to produce 



(2.27) Elog(l 



bL%) 



log (l + bu 2 
cosh 2 TTU 



du 



~(h(b,n)-4h(b,27r)) 
logi=-21og_L-^ 



1 

7h 



2 

71 



The duplication formula 

(2.28) ^(2z) = l 1 p(z) + ^(z + l)+log2, 

that appears as entry 8.365.6 in 7 , reduces (|2.27[) to the stated form. 



□ 
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3. The periodicity of the modified Bernoulli numbers Bfn+i 

This section uses the umbral method to express the generating function of the 
modified Bernoulli numbers B* t in terms of the digamma function ip(x). The period- 
icity of £>2n+i in Theorem 11.11 follows from this computation. Zagier [21] establishes 
this result by using the expression 

oo oo D r 

(3-1) 2 £ = J2 VTT^W ~ 2 l0g ^ " 

In the proofs given here the generating function employed admits an explicit ex- 
pression. 

Theorem 3.1. The generating function of the sequence {B^} is given by 

oo 

(3.2) F B , (z) := B* n z n = -- log z - (z + l/z - 1) . 



Proof. Start with 



OO OO 



n— 1 r— 

oo oo 



2r / n + r 



The second term is — log(l — z). Interchanging the order of summation in the first 
term gives 



oo oo 



n—1 r— 1 



n + r\ B r B r ^ n ( n + r-l)l 



2; 



2r n + r ^ (2r)\ ^ (n-r)\ 

/ r— 1 n— r 

and the inner sum is identified as 

" „(„ + r-l)I ~ „ +r (2r + m-l)!) ... : 



' ^ (n — r* 1 1 ^ — J m.i / 



(n-r)! m! '"' ^'(l- z ) 2r ' 

Therefore 

" B r 



(3.3) F B , (,) = - log(l - + £ ^ ' 



The rules of umbral calculus now give an expression for Fb* {z). The identity 

E— — =- = — eval < - log [ 1 



^2r(l_ z) - [2 ^ (l- Z ) 

gives 

(3.4) F B * 0) = -eval 1 1 log ( ( I - : r : 3? 



2 



s 
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Further reduction yields 

log ^(1 — z) 2 — 03z^ = logz + log 

= log z + log 

= log z + log 



(1-z) 2 



-03 



(i- z y 



03 + 1 



03 1 



(1-zY 



using (|2.9[) . The result now follows from Theorem 12.51 □ 

The generating function of -B^n+i is now obtained from Theorem 13. II The proof 
presented here is similar to the one given in [B] . 

Theorem 3.2. The generating function of the sequence of odd-order modified 
Bernoulli numbers is given by 

°° 9 7,5|3 q 

(3.5) Gb* (z) = £ B* n+1 z 2 ^ = 3Z Z 4 "/ 12 + _;, + Z ' 3Z . 
Proof. Start with 

(3.6) FB ' iz) ^ FB ' hz) =±B L+1 z^. 

To evaluate Fb*(—z) use the relation (|3.4I) and the operational rule from Lemma 
12.11 to obtain 

^2 



2F B . (-z) = -eval {log ((l + z) 2 + z03)} 



log z - eval i log ( - + ^ + 03 



log z — eval < log 



OS 



log z — tp 



Therefore 



G B * (z) 



F b * (z) - F B . {-z) 



"7 IV- (l + 



= iv(« + iA + 2)-^(« + i/«-i). 

Now use the relation (entry 8.365 in [7]) 

m — 1 

(3.7) tf(z + m) = + V — 

* — » •/ 4- 



fe=0 



z + fc 



to obtain the result. 



□ 
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Theorem 11.11 is now obtained as a consequence of Theorem 13.2 



Corollary 3.3. The sequence of odd-order modified Bernoulli numbers -B^n+i * s 
periodic of period 6. 

4. An inversion formula for the modified Bernoulli numbers 

This section discusses an expression for the classical Bernoulli numbers B n in 
terms of the modified ones S*. The result appears already in [21], but the proof 
presented here extends directly to the polynomial case as stated in Theorem 18.11 
The details are simplified by introducing a minor adjustment of B*. 

Lemma 4.1. Define B n = B* — 1/n. Then 



(4.1) 



5n = E 



fe=l 



n + k - 1\ B h 



- k 



2k 



Proof. The definition of B* in (| 1 .4[) produces 



k=l 



n + k\ B 



Then use 



1 



2k Jn + k 
(n + k-l)\ 



2k Jn + k (2k — l)!(n - k)\ 2k 



to deduce the claim. 

The inversion result is stated next. 



□ 



Theorem 4.2. The sequence of Bernoulli numbers B n are given in terms of the 
modified Bernoulli numbers B* by 



B n = 2n^(-1 



n+k 



k=l 



2n-l 

n — k 



2n- 1 
n — k — 1 



Bl 



•2(-iy 



2n- 1 
n 



for n > 1 . 

Proof. The inversion formulas 

f n + p + k 
n — k 



k=0 



b k , and b n = ^(-1) 



k+n 



k=0 



2n + p 
n — k 



2n + p 
n — k — 1 



ak, 



are given in [HI (23), p. 67]. Applying it to the sequence B n gives 



(4.2) 



fe = l 

The result now follows from 

~ '2n - 1 
n — /c 



(4-3) ^(-1 
fc=i 

To prove the identity 

1 

k 



2n- 1 

n — k 



2n- 1 

n — fc — 1 



2n - 1 

n — /c — 1 



wr 



ite the summand as 



(4.4) 



2n- 1 

n — k 



2n - 1 

n — A; — 1 



(-1)™ f2n-l 

7) 



2n 
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and convert (|4.3|) to 

£<-<-*) ~(V 

This follows directly from the basic sum 

□ 

5. A GENERATING FUNCTION FOR ZAGIER POLYNOMIALS 

This section gives the generating function of the Zagier polynomials 

n + r\ B r (x) 



2r I n + r 



(5-1) B*(x)=J2 

r=0 

The proof is similar to Theorem 13. 1[ so just an outline of the proof is presented. 
Theorem 5.1. The generating function of the sequence {B*(x)} is given by 

oo 

(5.2) F B . (x; z) = Y J B* n {x)z n = -- logz - (z + l/z-l-x). 

n=l 

Proof. The starting point is the polynomial variation of (|3.4j) in the form 



Fb* {x; z) — — eval 



(ibg((l-z) 2 -«»(»)) 



I 2 



= -eval | i log ((1 - zf -zx- z03) j 

= -eval | i log (l - 2z + z 2 - zx - zOS) 

= — — logz — eval / — log (l/z — 2 + z — x — 03) 
Now use — 03 = 03 + 1 to obtain 

(5.3) F B .(ir;z) = -ilogz- ieval{log (l/z + z - 1 - x + 03)} . 



The final claim now follows from Theorem 12.51 □ 
Corollary 5.2. The generating function of the sequence {(— l) n B* (x)} is given by 

oo 

(5.4) F B , (x; -z) := ^(-l) n B^(x)^ = -- log* - (z + l/z + 2 + *) . 

n=l 

Proof. Replacing z by — z in the third line of the proof of Theorem 15.11 gives 
F B »(x;-z) = -eval ji log (l + 2z + z 2 + zx + z03) 

= - i log z - eval | i log (z + l/z + 2 + a; + 03) | . 
As before, the result now comes from Theorem 12.51 □ 
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The next step is to provide analytic expressions for the generating functions 
of the subsequences {B% n+1 (x)} an d {^2n( a; )}- These formulas will be used in 
Section [T^] to obtain information about these subsequences and in particular to 
discuss periodic subsequences in Theorem II .21 



Corollary 5.3. The generating functions of the odd and even parts of the sequence 
of Zagier polynomials are given by 

oo 

(5.5) Yl B 2n+i(x)z 2n+1 = -^>{z + l/z + 2 + x)- -ip (z + l/z-l-x), 

n=0 



°° 1 1 1 

(5.6) J2 B 2n(x)z 2n = --logz--i,(z + l/z + 2 + x)- (z + l/z-l-x). 



The results of Corollarv l5.3l correspond to the analogue of the ordinary generating 
function for the Bernoulli polynomials. This is expressed in terms on the Hurwitz 
zeta function as stated in Theorem 15.41 The latter is defined by 

oo 

and it has the integral representation (see |18j . page 76) 

(5.8) ((s, a) = — - / —dt, Res>l, Rea>0. 

r(s) Jo 1 -e 

The exponential generating function (jl.6l) . because it is given by an elementary 
function, is employed more frequently. 

Theorem 5.4. The generating function of the Bernoulli polynomials B n {x) is given 
by 

OO ^ 

(5.9) YB n (x)z n = -C(2,l/z-x + l). 
t-^ z 

Proof. The integral representation of the gamma function 

/•OO 

(5.10) r(s) = / u^e^du, Res>0 

Jo 

and the special value T(n + 1) = nl give 

00 POO 00 R ( \ 

(5.11) YB n {x)z n = e- u Y^^(zu) n du. 

rr£ Jo n - 



n=0 ° u n=0 

The generating function ()1.6j) is used to produce 
(5.12) 5>„(*),» = 2 / _ — 

m — n ^ 



n=0 

The change of variables v — zu and (15.81) complete the proof. □ 
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6. Some arithmetic questions 

There is marked difference in the arithmetical behavior of the numbers B* L (j) 
according to the parity of n. For instance 



{B*(0) : 1 < n < 10} = 
and 



3 1 1 27 1 29 1 451 1 65 1 

4' 24' _ 4' ^80' ~4' ~1260' 4' 1120 I' ~ 264 J 



, 5 25 5 133 9 3751 15 4931 19 12551 
{B*(l): l< n <W}={-, 



4' 24' 4' 80 ' 4' 1260' 4 ' 1120 4 ' 264 J ' 
On the other hand, keeping n fixed and varying j gives 



{bUj)-- i<j<io} 

and 



5 7 9 11 13 15 17 19 21 23 

I' I' I' T' T' T' T'TI'I 



is*n i< • < ioi - 1 25 61 109 169 241 325 421 529 649 781 \ 

{ 2 {3) ■ <J< l ~ I 24' 24' -24> -24~> "^f "24" "24"' "24J ' 

This suggests that every element in the list {B*(j) : j > 1} has a denominator 
that is independent of j, therefore this value is also the denominator of the modified 
Bernoulli number B*. Assume that this is true and define a(n) be this common 
value; that is, 

(6.1) a(n) = denominator^*). 

As usual, the parity of n plays a role in the results. 

The next theorem shows, for the case n odd, that the function a(n) is well 
defined. 

Theorem 6.1. For j S Z, the values 4i?2„+i(j) ar £ integers. That is, 

(6.2) a(2n + l)=4. 
Proof. The generating function (|5 .3[) gives 



J2 4S 2 *„+iO> 2n+1 = ip(z+ 1/z + j + 2) ~ tp (z + 1/z 

n=0 

= i>(z + l/z) + j: 1 -iP(z + l/z) + 

k=0 1 
J 1 



k=o z + l/z-j-l + k 

J+l 1 3+1 1 



z + l/z + fc ^z + l/z-j-l + fc z + l/z 
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Replace k by j + 1 — k in the second sum to obtain 

1 



5>B 2 *„ + iO> 2n+1 



E 



fc=0 

J'+l 

2 *E 

fc=0 

2 -E 



z + l/z + fc z + l/z 



z 2 + 1 



(z 2 + l) 2 - k 2 Z 2 
z 2 + 1 



(Z 2 + l) 2 - k 2 Z 2 



Z 2 + l 



Z 2 + l 



This implies 4 J B*„ +1 (j) e Z. 



z + l/z 



□ 



Note 6.2. Arithmetic questions for the numbers B| n (j) seem to be more delicate. 
The values a(2n) seem to be divisible by 4 and the list of ja(2n) begins with 

{6, 20, 315, 280, 66, 3003, 78, 9520, 305235, 20900, 138, 19734, 6, 7540, 15575175}, 

for 1 < n < 15. The exact power of a prime p that divides a(2n) exhibits some 
interesting patterns. For instance, Figure [T] shows this function for p = 2. 





20 40 



100 120 



Figure 1. Power of 2 that divides denominator of S|„(j) 



The data suggests that the prime factors of a(2n) are bounded by 2n + 1. These 
questions will be addressed in a future paper. 



7. Some auxiliary binomial sums 

This section contains the proofs of two identities for some sums involving bino- 
mial coefficients. These sums will be used in the Section [5] to give an expression 
for the derivatives of Zagier polynomials. The identities given here are established 
using the method of creative telescoping described in |13j . 



Lemma 7.1. For nGN 

(7.i) ]T(-iy 



r 2(r + 1) fir - 1\ (n + r + 1 
n + r + ll r )\2r + 2 
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Proof. The summand on the left-hand side is written as 

(-l) r (n + r)\ 



F(n,r) 



2{2r + l)rl 2 (n-r-l)\' 



Observe that F(n, r) vanishes for r < or r > n— 1. The method of Wilf-Zeilberger 
lends the companion function 

(7.2) G(n,r)- 



(n + l)(r - l)! 2 (n - r + 1)! 
together with the second order recurrence 

(7.3) F(n + 2, r) - F{n, r) = G(n, r + 1) - G(n, r). 

Sum both sides over all integers r and check that the right-hand sum vanishes to 
produce 

(7.4) ^F(n + 2,r) =^F(n,r). 
Dehne 

n-l 

(7.5) f(n) = ^F(n,r). 

r=0 

Then CLl]) gives /(n + 2) = f{n). The initial conditions /(l) = 1/2 and /(2) = 
show that 



(7.6) f(n) 



2 

for n even. 



., for n odd 



The desired sum starts at r = 1, so its value is /(n) — -F(n, 0). Thus, F(n, 0) = n/2 
gives the result. □ 

Lemma 7.2. For neN and 1 < k < n — 1, define 

lnn , , ^2(-l) r r(r + l)/n + r + l\ T/2r-l\ / 2r - 1 

(7.7) M (n,fc) = JJ 



n + r + 1 V 2r + 2 

-k 



r ■ 



— k J \r — k — 1 



T/ien for n even 

(7.8) w(n,fc) = < 
and for n odd 

(7.9) u(n,k) = 



—k for k odd, 
for k even 



for k odd, 
k for k even. 



Proof. A routine binomial simplification gives 

n-l 

(7.10) u(n,k) = fc^(-l) r 



n-l / | \ / o 

n + r \ / Zr 



r—k 



2r + l\r-k 
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This motivates the definition u(n,k) = u(n,k)/k and the assertion in (|7.8j) and 
(|7.9|) amounts to showing 



!+l for n odd, k even 
— 1 for n even, k odd 
otherwise. 

The proof is similar to the one presented for Lemma 17.11 Introduce the functions 
F(n,r,k) = (-l) r ( 2 "+ r i) ( r 2 - k ) and use the WZ-method to find the function 

(7-12) G (n ) r )fc ) = t + , 1)(2 V 1)( : 1) Tf n + r> )f V 

v ' y ' ' ' (n + k + l)(n-k + l) \2r - lj \r - k - 1, 

companion to F and the equation 

(7.13) F(n + 2, r, k) - F(n, r, k) = G(n, r + 1, k) - G(n, r, k). 

The argument is completed as before. □ 

8. The derivatives of Zagier polynomials 

Differentiation of the generating function for Bernoulli polynomials (|1.6j) gives 
the relation 

(8.1) — B n {x) =nfl n _i(x). 

This section presents the analogous result for the Zagier polynomials. The proof 
employs an expression for B n (x) in terms of £?*(x); that is the inversion of (jl.5l) . 
The proof is identical to that of Theorem 14.21 so it is omitted. 

Theorem 8.1. The sequence of Bernoulli polynomial B n (x) is given in terms of 
the Zagier polynomials B* t (x) by 



B n (x) = 2nJ2(-l) 



2n - 1\ / 2n - 1 

n — k I \n — k — 1 



B* k (x) + 2(-l) n 



2n- 1 

n 



fe=l 

for n > 1. 

The analogue of (|8.ip is established next. 
Theorem 8.2. TTie derivatives of the Zagier polynomials satisfy the relation 

LfJ 

R*W = 

dx 



^B* n (x) = £(2j ~ ^)B* 2j -x(x) for n even 



and 

LfJ 



d_ 

~dx~ ny ~' 2 



B*(x) = - + 2 J B 2j( X ) f° r n odd - 



3=1 

Proof. Differentiating (|5.1[) and using (|8.1[) gives 

nV y ^ V 2r + 2 Jn + r + l y ' 

r— x ' 

= 2+M 2r + 2 J^+rTT ^ 

r— 1 x / 



1(5 



A. DIXIT, V. H. MOLL, AND C. VIGNAT 



for n > 1. The sum above (without the term n/2) is transformed using Theorem 
18.11 to produce 



1 r 2r — 1\ /n + r + 1 



r=0 
n-1 



2r + 2 /n + r + 1 



Kr + l) /n + r + 1\ ^ ( _ 1)k +r 
^n + r + ll 2r + 2 / ^ l ' 

r=l v 7 fc=l 



2r - 1\ / 2r - 1 

r — k 1 \r — fc — 1 



Denote the first sum by Si and the second one by S^- 

Lemma f7.ll shows that Si = — I §J . To evaluate S 2 , reverse the order of sum- 
mation to obtain 

n-i 

(8.3) S 2 = £(-l) fc u (n,*0^(z), 

k=i 

with u(n,k) defined in (|7.7j) . This gives 

, n— 1 

(8-4) -B„(z) = \ [l\ +^(-l) fe U (n,fc)^(x), 

and the proof now follows from the values of u(n, k) given in Lemma 17.21 □ 

9. Some basics on Chebyshev polynomials 

This section contains some basic information about the Chebyshev polynomials 
of first and second kind, denoted by T n (x) and U n (x), respectively. These properties 
will be used to establish some results on Zagier polynomials and the relation between 
these two families of polynomials will be clarified in Section [TOl 

The Chebsyhev polynomials of the first kind T„ are defined by 

(9.1) T n (cos 9) = cos n9, n E N U {0} 

and the companion Chebyshev polynomial of the second kind U n by 

sin((n + 1)9) 

9.2 U n cos 9) = VV . ' ' . 

smfci 

Their generating functions are given by 

oo 

(9.3) Y, T ^ tn 

n=0 

and 

oo 

(9.4) ^f/„(i)i" 
The even part of this series, given by 

°° 1 4- t 2 

(9-5) E^™(^ 2n = IT^ 



l-xt 



l-2xt + t 2 



1 



1 - 2xt + t 2 ' 



. _, . 2x 2 )t 2 +t 4 '' 

n— v ' 

will be used in Section [TTJ Many properties of these families may be found in |12) . 
The formulas for the generating functions also appear in [16l 22 : 3 : 8, page 199]. 
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Differentiation of the generating function for T n (x) gives the basic identity 



(9.6) 

The expression 
(9.7) 



—T n {x) = nU n -i(x). 
ax 



U n (x) 



(x + Vx 2 - l) n+1 - (x - Vx 2 - l) n+1 



2Vx 2 ^l 

will be used in the arguments presented below. 

The Chebyshev polynomials have a hypergeometric representation in the form 

(9.8) T n ( X ) = 2 F l (_„,„; 

U n (x) = ( n +l) 2 F 1 (-n,n+2;§;^). 

These appear in jTSJ p. 394, equation (15.9.5) and(15.9.6)]. 

The next statement is an expression of the Chebyshev polynomial T n that will be 
used to establish, in Theorem lll.il a symmetry property of the Zagier polynomials. 

Lemma 9.1. The Chebyshev polynomial T n {x) satisfies 

'n + r\ x r 



(9.9) 



E 

!-=0 



2r J n + r n 



^(1 + 



T » 1.2 + 1 



Proof. The representation (JH]) and (|) r = 2~ 2r (2r)!/r! yield 
= 2-Fi {-n, n; |, -f ) 

r=0 \1> r 
n 

= ^2{n(n - 1) ■■•(»- (r - 1))} {n(n + 1) • • • (n + r - 1)} 



r=0 



(2r)\ 



= »E 

^ (n + r-1)! r 
j^(2r)l(n-r)! 

This verifies the claim. 



\ (n - r)\ {(n - r + l)(n - r + 2) • • • (n + r - 1)} 



(2r)!(n-r)! 



□ 



Lemma 9.2. The Zagier polynomial B^(x) is related to the Chebyshev polynomial 
T n (x) via 



Proof. This is simply the umbral version of (19.91) . 
Lemma 9.3. The Zagier polynomial B*(x) is given by 



□ 



(9.10) 



B*(x) = — E 



. x 1 

2 + 2 iLi 



Proof. The result now follows from Lemma 19.21 Theorem 12.31 and the umbral rule 
(1231). □ 
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The result of Lemma 19.31 is now used to give an umbral proof of Theorem 18.21 

Proof. The computation is simpler with B n (x) — B* n (x — |). Differentiate the 
statement of Lemma 19.31 to obtain 



(9.11) 



^B n (x) = — E 
ax 2n 



tL 



x + iLb 



1. 



-E 



In the case of even degree, this gives 



(9.12) 

and using the identity 
(9.13) 



d - , N 1 

-rB 2n {x) = -E 

ax 2 



U 2n -i 



2 

x + iLb 



x + iLf 



U 2n -i(x) = 2y^T 2fc -i(a 



fc=i 



that is entry 18.18.33 in [12] . it follows that 

x + iLb 



(9.14) ^-B 2n (x) 
ax 



zJ T2k - 1 

Lfc=l 



k=l 



as claimed. The same argument works for n odd using 12, 18.18.32]: 

n 

(9.15) U 2n (x) =2j2T2 k {x)- 1. 



fe=0 



In [: 
(10.1) 



10. The Zagier-Chebyshev connection 
after the proof of the identity 



r=0 



n + r\ B 



>2r 



1r J n + r ' 



□ 



the author remarks that the second term has a pleasing similarity to the equation 
(|1.4[) . This section contains a representation of the Zagier polynomials B^(x) in 
terms of the Chebyshev polynomials of the second kind U n (x). The expressions con- 
tain terms that also have pleasing similarity to the definition of Zagier polynomials. 
The results are naturally divided according to the parity of n. 



Theorem 10.1. The Zagier polynomials are given by 



(10.2) 2BUx) = 



n+r 



r=0 



n + r\ B 2r (x) 



2r I n + r 



U 2n -\ 



(I) 



+ U 2n . 



x + 1 



(10.3) 2BZ n+1 (x)=J2(- i y l+r 



n + r + l\B 



>2r+l 



(X) 



2r + l n + r + 1 



-u 2 . 



(!)■ 



■U 2 



X+1 



Proof. The proof is presented for the even degree case, the proof is similar for odd 
degree. 

Theorem 15.11 gives the generating function for B^(x). Its even part yields 



2j2 B L(x)z 2n = --logz--i ) (l/z+z-x-l)--log{-z)--i>(-l/z-z-x-l). 
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The functional equation ip(t + 1) = ijj{t) +l/t gives 

oo 

(10.4) 2 Y, B* 2n {x)z 2n = H(x, z) + H{x, -z) 



n=l 
1 



2 \l/z + z + x + 3 l/z + z + x + 2 l/z + z-x-3 l/z + z-x-2 

with 

(10.5) H(x, z) = -- (logz + ^(l/z + z - x - 3)) . 



The umbral method and Theorem 12 . 5 1 give 

2H(x,z) = -logz-eval(log(l/z + z-x-3 + < B)) 
= — log z — eval (log(l/z + z — x — 4 — 25)) 
= —eval (log(l + z 2 — 4z — zx — cc23)) 

= -log(l + z 2 )-eval(log(l- Z ^ + 4) )) 

, ,1 2n , /v 2 - (z $ B(a; + 4))'-\ 
= - log(l + z 2 ) + eval 2^ v ' 



log(l + z 2 )+]T 



\r=l 

5,.(x + 4)z 



r(l + z 2 ) r 



r (l +z 2)r 
r— 1 v ' 

Therefore 

H(x, z) + H(x, -z) = log(l + z 2 )+Y 2 r(l + z 2 V ' 

r— 1 ^ ' 

Now observe that 

y, £ 2r Qr + 4)z 2 '- = y, z 2 '-B 2r (x + 4) y, (2r) n (- z 2 ) 
^ 2r(l + z 2 ) 2 '' ^ 2r ^ n\ 

r—1 x ' 



n=0 

(-l) r B 2 Ax + 4) yv_ iyi (n + r-l 
^ 2r ^ [ ' V 2r-l 

r—1 n— r v 

oo / n 



z 2n 



= E(E(-i) 

This produces 



/n + r\ S 2r (.T + 4) 

n— 1 \r— 1 



2r / n + r 



z 2 ™. 



The equation (|10.4[) now gives 



oo I n 



(10.6) 2YBL(^ 2n = Y[E(- i r +rfn + ^ B2AX + A) 



2r I n + r 



z 2n 



n—1 ri— 1 \r— 

l/i l l i N 

2 \l/z + z + i + 3 l/z + z + a; + 2 l/z + z-x-3 l/z + z-x-2, 
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The generating function for U n -i(x), given in (|9.4|) . is written as 

oo 1 

and the rational function appearing in (|10.6[) is expressed as 
^ / f-x-3\ TT f- x -2\ TT fx + 3\ TT fx + 2 

£ K-iHr- -^Hr- Hn 



n=l v v ✓ v / v ' 

Using the fact that U n {x) has the same parity as n it simplifies to 



v — / / —x — 3\ /—a; — 2 



z 2 ". 



Therefore 

0R , n v\ n n+r^ + A -B2r(a: + 4) f -x-3 ^ f -x-2 

2B ^)=}J-V [ 2r J n + r + t/2 "- 1 ( v ^r-J +t/2 "- 1 ( v ^- 

Finally, replace x by — x — 3, use Theorem 111.11 and the symmetry of Bernoulli 
polynomials i?2n(l — x) = B2n(x), to obtain the result. □ 

The next result gives a representation for the difference of two Zagier polynomials 
in terms of the Chebyshev polynomials U n (x). 

Lemma 10.2. The Zagier polynomials satisfy 

(10.8) B* n (x + 1) = B*(x) + hj n - x (| + l) . 
It can be extended to 

(10.9) B*(x) - B*(x -k) = \jz Un-x + 1 

3=1 V 

Proof. Apply Lemma [2~2l and the representation of B^(x) in Lemma [9721 □ 
Note 10.3. The sum in (IT0~2)) equals 

(10.10) X>l)" +r (\ + ^ = 2B *^ X - 2 ^ 

\ 2r / n + r 

r=0 V 7 

The proof of this fact is given below. First, use it to express (I10.2[) as 

(10.11) B* 2n (x) B* n {x - 2) = 1 (u^-i (f) + Cbn-i ( ' + ' 

In this form, it can be extended directly to 
(10.12) 

B* 2n (x) - B* 2n (x -2k) = \ ^»-t (^) + U^-i ( X ~ 2j + 1 

The proof of (jlO.101) is given next. The identity 
(10.13) T 2n (x) = (-l) n T„(l - 2x 2 ) 
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that appears in [4] 7.2.10(7), page 550] is used in the proof. Start with 



r=0 \ / ' r=0 \ / 



(-» 2 (x))' 



n + r 



eval 



eval 



(-1)" ( <&\x) 

~ -L i 



n 

(-i)" 

n 



T n 1 - 2 



2eVal ^ 2^ T2 " 



2eval<! — T 2n 



1 /»(*)- 2 



2eval<i — T 2n 



1 /®(i-2) 



2n + r\ 2) 
2r / 2n + r 



- 2S 2 *„(x-2). 

A special case of Theorem 110.11 gives a simple proof of (|10.1[) . 
Corollary 10.4. The modified Bernoulli numbers B^ are given by 



(10.14) 



2B% 



_ fn + r\ B 



r=0 



2r I n + r 



Here (^) is £/ie Jacobi symbol. 

Proof. Let x = in Theorem ll0.ll use the value U2 n -i(0) = and observe that 

1 if n = 1 mod 3, 



(10.15) 



C/ 2n -i (|) = 



— 1 if n — —1 mod 3, 
if n = mod 3, 



can be written as U^n-i (5) — ( ~S ) 

The next result appears in |21j . 
Corollary 10.5. Let n e N. TAen 



□ 



(10.16) 



Bl 



In 

E 

r=0 



2n + r\ {-l) r B r 



2r J 2n + r 



Proof. The right-hand side of (|10.16[) is S| n (l). Therefore, the statement becomes 
B*„ + n = i?2„(l)- This is established by letting x = in l|10.8p and the value 



J 2n -r n 
(10.17) 



. . sin2n0 
(V2n-i(l) = hm — : — — = 2n. 



>o sm I 



□ 



Theorem 110.11 is now used to produce another proof of Theorem 11.11 
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Corollary 10.6. The modified Bernoulli numbers B 2n , j are given by 
(10.18) B 2n+1 - + -U 2n { 2 ) - + sin I I . 

In particular, B 2n+1 is periodic of period 6. 

Proof. Put x = in (|10.3[) and observe that only the term r — survives in the sum. 
Now use the value U 2n (0) = (-1)" and let 9 = tt/3 in (J9T2J) to get the result. □ 

Corollary 10.7. For n € N 

(10-19) 2B 2 *„ +1 (|) = tf 2n (i) + C/ 2n (|) . 

Proo/. Let x = | in (|TUT5j) and use the fact that Bj{\) = -(1 - 2 1 ~ j )B j . □ 

11. A REFLECTION SYMMETRY OF THE ZAGIER POLYNOMIALS 

The classical Bernoulli polynomials B n (x) exhibit symmetry with respect to the 
line x = i in the form 

(11.1) B n (l-x) = (-l) n B n (x). 

This section describes the corresponding property for the Zagier polynomials: their 
symmetry is with respect to the line x = — %. 

Theorem 11.1. The Zagier polynomials satisfy the relation 

(11.2) B*(-x-3) = (-l) n B*(x). 

Proof. The first proof uses the generating function Fb*(x,z). Replacing (x, z) by 
(—a; — 3, — z) in the second line of the proof of Lemma 15 . 1 1 gives 



Fg* (—x — 3, —z) = -cval|ilog((l + z) 2 -z(x + 3) + z < B) 

= — — log z — eval | — log (z+l/z + *8 — 1 — x) 
= F B *(x,z). 

This proves the statement. □ 

A second proof of Theorem 111.11 uses the expression for the Zagier polynomi- 
als B^(x) in terms of the Chebyshev polynomials T n (x). Indeed, using T n (z) — 

(-l) n Tn(z), 

B n(-x ~ 3 = -T n + 1 

n \ 2 

1 f-x + K 1 

J. 7 



n \ 2 2 

(-1)™ fx-<s 

n 

(-1)" (x + % + 1 ( 1 



2 


r 2 


x + < B H 


- 1 


2 






f 1 



n \ 2 2 

(-1)" T 

J- n. 



-l) n B* n {x). 
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The rest of the section gives a third proof of Theorem lll.il 

Proof. The induction hypothesis states that B^(—x — 3) = (-l) m B* n (x) for all 
m < n. The discussion is divided according to the parity of n. 
Case 1. For n is even, Theorem 18.21 gives 

d n/2 

—B* n (-x-3) = -^(2j-l) J B| i _ 1 (-x-3) 

3=1 

n/2 

= -^(2j-l)(-l) 2 i- 1 B| J ._ 1 (a ; ) 

3=1 
n/2 

3 = 1 

ax 

It follows that -B„(— x — 3) and B^(x) differ by a constant. Now evaluate at a? = — | 
to see that this constant vanishes. 

Case 2. Now assume n is odd. The previous argument now gives 

(11.3) B* n (-x-3) = -B* n (x)+C n 

for some constant C n . It remains to show C n = 0. 
Iterating the relation 

(11.4) B n (x + 1) = B n {x)+nx n - 1 
gives 

(11.5) B n (x + 3) = B n (x) + nx n - 1 + n(x + l)™" 1 + n(x + 2)"" 1 . 
Replace x = — | in (lll.3[) and in (|1 1 .5[) and observe that 

(H-6) C» = 2 E 7TT7 * (|) - r - r (- J)"" 

r=0 

Thus, to show C n = 0, it is required to prove 

n fn+r\ n fn+rS 

The left-hand side is nothing but B* (|). The right-hand side is K(-|)+V^(-|), 
where 



(11.8) K(x) = E 



n fn+r* 
2r , 



n + r 

r=0 



Lemma EH] shows that 

(11.9) V n (x) = ~T n (~ + l). 

Hence it suffices to show that 

(11-10) 2B* (|) = [/„_! (i) + l/ n _i (|) . 

This is the result of Corollary [11171 □ 
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Note 11.2. We note that unlike (fTCL2]) . the formula in (fTTH]) . of which (jll.lOjl is a 

special case, does not use the symmetry B 2n+1 {— x — 3) = — B 2n+1 (x) in its proof. 

12. Values of Zagier polynomials that yield periodic sequences 

The original observation of Zagier, that B 2n+1 — -^2n+i(^) is a periodic sequence 
(with period 6 and values {§, —i, — i, i, i, — |}) is extended here to other values 
of x. The first part of the discussion is to show that periodicity of B 2n+1 (x) implies 
that 2x is an integer. 

The discussion begins with an elementary statement. 

Lemma 12.1. The sequence {a n } is periodic, with minimal period p, if and only 
if its generating function 



oo 

(12.1) A(z) = J2 a " zT ' 

n=0 



is a rational function of z such that, when written in reduced form, the denominator 
has the form D(z) — 1 — z p . 

Special values of S| n+1 (cc). The case considered here discusses values of x such 
that {B 2n+1 (x)} i s a periodic sequence. The generating function of this sequence 
is given in (|5.5[) as 



1 1 

(12.2) B* 2n+l (x)z 2n+1 = (z + l/z + 2 + x)- -V (z + l/z-l-x). 

n=0 

Proposition 12.2. Let b e R be fixed. Then 

(12.3) ip(t + b)-ip(t) = R(t) 
for some rational function R(t) if and only if b € Z. 

Proof. Assume b € Z. It is clear that b may be assumed positive. Iteration of 
*/>(*+ 1) = tp(t) + l/t yields 

6-1 

(12.4) ip(t + b)=m + J2 

k=0 

Therefore ip(t + b) — ip(t) is a rational function. To prove the converse, assume 
(|12.3|) holds for some rational function R. Integrating both sides with respect to t 
gives 

(12.5) hxT(t + b) - lnr(t) = R^t) + lnR 2 (t) + C x 

for a pair of rational functions R\, Ri (coming from the integration of R(t)) and a 
constant of integration C\ . It follows that 

(12 6) F( * + b) - e R ^ 

[ ' C 2 R 2 (t)T(t)- e ■ 

The singularities of the left-hand side are (at most) poles. On the other hand, the 
presence of a pole of i?i (t) produces an essential singularity for the right-hand side 
of (|12.6I) . It follows that Ri(t) is a polynomial. Comparing the behavior of (|12.6[) 
as t — > ±oo shows that R± must be a constant; that is, 

(12.7) T{t + b) = C 3 R 2 (t)T(t). 
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The set equality 

(12.8) {b- k : k g N} = {-k : k g N}U{t x , i 2 , • • • ,t r } 

where t% are the poles of R comes from comparing poles in (| 12 . T|) . Now take k g N 
sufhciently large so that b — k ^ ti. Then b — k = —k\ for some &i g N. Therefore 
b = k — fci G Z, as claimed. □ 

The next lemma deals with the transition from the variable z to z + 1/z. 

Lemma 12.3. Assume R(z) is a rational function that satisfies R(z) — R(l/z). 
Then R is a function of 1/z + z only. 

Proof. It is assumed that 

(12 9) R(z) - a ° + aiZ H a " Z " - a o + n i/ z j !: a "/ z " 

bo + M + • • • + bm-z m ~~ b Q + bi/z + ■ ■ ■ + b m /z m ' 

Now use the fact that 

,, 01 m " ^ • r u U u+U 

(12.10) — = — implies — = — = 

v V v V v + V 

to conclude that 

2a + oi (z + 1/z) + a 2 (z 2 + 1/z 2 ) + ■■■ + a n (z n + l/z n ) 



R(z) 



2b + h(z+ 1/z) + b 2 (z 2 + 1/z 2 ) + ■■■ + b rn (z m + l/z m ) ' 

The conclusion follows from the fact that z J + 1/z 3 is a polynomial in z + l/z. This 
is given in entry 1.331.3 of [7]. □ 

Theorem 12.4. The generating function 

oo 

(12.11) E^+iW z2n+1 

is a rational function of z if and only if 2x G Z. 

Proof. Assume (|12.11l) is a rational function of z. Then (|5.5[) implies that 

(12.12) Tj}(z + l/z + 2 + x)-ip(z + l/z-l- x) = A(z) 

with A a rational function of z. The left-hand side of (|12.12[) is invariant under 
z i — ^ 1/z, therefore Lemma [12.31 shows that A(z) = _B(z + 1/z), for some rational 
function B. Now rewrite (|12.12[) as 

(12.13) ip(t + 2x + 3) - i/}(t) = B(t + 1 + x) 

with £ = z + l/z — 1 — x. Proposition 112. 21 shows that 2x G Z. 

To establish the converse, assume 2x G Z. The identity (|5.5p shows that 

oo 

(12.14) 4 B* 2n+1 (x)z 2n+1 = i/>Qt + 2x + 3) - 4>{t) 

ra=0 

with i = z + 1/z — 1 — x. Proposition 112.21 shows that + 2x + 3) — ip(t) is a 
rational function of £ and hence a rational function of z. □ 

Corollary 12.5. Assume the sequence {-B2„ +1 (x)} zs periodic. Then 2x G Z. 

Proof. The hypothesis implies that the generating function in (|12.11j) is a rational 
function. Theorem 112.41 gives the conclusion. □ 
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The quest for values of x that produce periodic sequences B 2n +i(x) l& now re ~ 
duced to the set Z U (Z + h). The symmetry given in Theorem 111.11 implies that 



one may assume x < 



2 ' 



12.1. Integer values of x. The nature of the sequence {B 2n+ i(x)} is discussed 
next for x = k £ Z. 



Theorem 12.6. Let n G N and fc > 3. Then 



1 1 fe ~ 2 / ' 

(12.15) S| n+1 (-fc) = --LT 2n (0) - - U *n [% 



Proof. This is just a special case of (| 10 .9[) with x = 0. Use (|10.18[) and the fact 
that U 2n (0) = (-1)™. □ 

The next step is to show that {-B 2n _|_i(— fc)} is not periodic for k > 5. 

Lemma 12.7. Assume j > 3. Then f7 2n (§) > 0. 

Proof. This comes directly from (|9.7I) . □ 
Proposition 12.8. The sequence {-B2 I j +1 (— fc)} is no£ periodic for k > 5. 
Proof. The identity (112.15)) is written as 

-2£ 2 * n+1 (-fc) = + U 2n (l/2) + U 2n (l) + U 2n (3/2) + J2 U 2n (f 

3 '=4 ^ 

> + C/ 2 „(l/2) + C/ 2 „(l) + C/ 2 „(3/2). 

The value 



(12.16) C/ 2 „ 



2/ 



/• ^\ 2n+l / _\ 2n+l 

3 + V5\ /3-V5\ 



shows that |-B^„ +1 (-fc)} is not bounded. To obtain (|12.16[) . use x = § in (j9?7Tl . □ 

The next result shows that, after a linear modification, the case k = —4 produces 
another periodic example. 

Proposition 12.9. The sequence {-B 2 „_|-i( — 4) + n} is 6-periodic. 
Proof. The value /c = 4 in (|12.15|) gives 

(12.17) S 2 * n+ i(-4) = -|c/ 2 „(0) - ^E/ 2 „ (|) - i[/ 2 „(l). 
The values C/ 2 n(0) = (—1)™ is 2-periodic and 

(12.18) [/ 2 „(|) = -^sin((2n + l)7r/3) 

is 3-periodic (with values 0, —1, +1). The expression (|9.2p . in the limit as 8 — > 0, 
gives C/2n(l) = 2n + 1. The proof is complete. □ 

Corollary 12.10. The sequence {-B2 ?l+1 (— 3)} is 6-periodic. 
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Proof. Choose k — 3 in (112. 15|) to obtain 

(12.19) B* 2n+1 (-3) = -\u 2n (0) - \u 2n Q 

As in Proposition 112.91 £/2n(0) is of period 2 and U-z n (|) is of period 3. □ 
Proposition 12.11. The sequence {-Bln+iC - 2)} is 2-periodic: 

(12-20) B« (_ 2= ^J . 



Proo/. Let fc = 2 in Theorem 1121)1 □ 

The rest of the integer values x are obtained by the symmetry rule given in 
Theorem 111.11 The study of the structure of the sequences F2„ +1 (fc) has been 
completed. The details are summarized in the next statement. 

Theorem 12.12. Let k € Z. Then 

a) {B2 n+1 (k)} is exponentially unbounded if k > 2 or k < —5; 

b) {i?2n+i(^) + n } * s ^-periodic for k = —4 or k = 1; 

c) {-E>2„ +1 (fc)} is ^-periodic if k = —3 or k = 0; 
efj {S2„ +1 (fc)} is 2-periodic if k — —2 or k = —1. 



12.2. Values of x G i + Z. The example x = — | is considered first. 
Proposition 12.13. For n e N U {0} ; 

3t 



(12-21) -| =0 



Proof. Theorem 111.11 states that F2 n+1 (— x — 3) = — B2 n+1 (x). Replacing x = — | 
gives the result. □ 

The symmetry of B2 n+1 (x) about x = — | shows that it suffices to consider 
values of the form k + \ for k > — 1. 

Theorem 12.14. For aZZ fc > -1, 

(12.22) BL+1 ( k+ l ) = lf:U 2 J?^±l). 

Proof. The proof is similar to that of Theorem ll2.6[ so the details are omitted. □ 
The next lemma produces an unbounded value in the sum (|12.22l) when k > 1. 
Lemma 12.15. For n e N 



22n+2 _ 2^2™ 



(12.23) U 2nKtil - ., 

Proof. This comes directly from (|9.7I) . □ 

The next examples deal with values of F2„ +1 (fc + |) that do not contain the 
unbounded term Uin (f)- 

Lemma 12.16. The sequence B'2 n+ i (— k) * s periodic. 
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Proof. Theorem [127141 with k = -1, and (|972|) give 

B* 2n+ i H) = l U *n (I) = -|= sin ((2n + 1)0) 



with cos^ = j. It follows from here that {B 2n+1 (— ^)} is not periodic. Indeed, if p 
were a period, then B% n+2p+1 (-§) = SJn+i (-|) implies 

(12.24) tan((2n + 1)0) = cotp0 for all n £ N. 

Thus 30 and must differ by an integer multiple of 7r; that is 20 = 7rm. This is 
impossible if cos0 = i. □ 

Lemma 12.17. The sequence Bfn+i (i) * s not periodic. 
Proof. In the case fc = 0, Theorem ll2.14l gives 

(12.25) B 2 *„ +1 (|) = \ [U 2n (|) + f/ 2 „ (I)] . 

To check that this is not a periodic sequence, use (|9.5I) to produce 

(12.26) ^^„U)+C/ 2 „(l)J* " 16*4 + 24*3 + 25*2 + 24* +16' 

n— 

Periodicity of B|n+i (^) implies that the poles of of the right-hand side in (|12.26[) 
must be roots of a polynomial of the form 1 — t p . In particular, the arguments of 
these poles must be rational multiples of it. One of these poles is t = (l + 3\/7i)/8, 
with argument a — cos -1 (~j. Therefore a must be a rational multiple of 7r. To 
obtain a contradiction, observe that 

/7rm\ 

(12.27) a; m) „ :=2cos(— j 

is a root of the monic polynomial 2T n (x/2). It follows that w m>n is an algebraic 
integer and a rational number (namely i). This implies that it must be an integer 
(see [T71 page 50]). This is a contradiction. □ 

These results are summarized in the next theorem. 

Theorem 12.18. There is no integer value of k ^ —2 for which {B 2n+1 (k + h)} 
is periodic. 

Special values of B 2n (x). The second case considered here deals with values 
of the subsequence B 2n {x). Symbolic experiments were unable to produce nice 
closed- forms for special values of B2 n (x), but the identity 

(12-28) B* 2n (-1) = B* 2n (-2) 

motivated the definition of the function 

(12.29) A* 2n (u) := £?*„(- 1 - u) - B* 2n (-l), for u£l. 

Lemma 12.19. For n £ N, the function A 2n (u) satisfies A 2n (u) — A 2n (l - u). 
Therefore, it suffices to describe A 2n (u) for u > 1. 

Proof. This follows directly from Theorem lll.il □ 

The next statement expresses the function A 2n in terms of the Chebsyshev poly- 
nomials of the second kind U2n-i{x). 
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Proposition 12.20. The function A\ n is given by 

, u+l 

(12.30) A* 2n {u) = -^2^1 

Proof. Iterate the identity (flTTS!) . □ 

The expression in (|5.6I) yields the next result. 

Lemma 12.21. The generating function of the sequence {£>2ri( x ) — B% n {— 1)} sat- 
isfies 

n—l 

with w — z + 1/z. 

The proof of the next result is similar to that of Theorem 112.41 
Corollary 12.22. The generating function 

oo 

(12-31) El^W-^J-l)]^ 

n=l 

is a rational function of z if and only if 2x G Z. 

The next statement is an analogue of Theorems 112.121 and 112.181 

Theorem 12.23. Let A* 2n {x) = B^J-l - x) - B^ n (-l) as above. Then 

1) The sequences A^^l) and A2 n (0) vanish identically. 

2) The sequences A2„(2) and At,^— 1) are periodic with period 3. The repeating 
values are {i, —5, 0}. 

3) The sequences ^4.2„(3) and ^2„(— 2) grow linearly in n. Moreover, ^2^(3) — n 
and A2 n (— 2) — n are periodic with period 3. The repeating values are {i, — ~, 0}. 

4) The sequence Ai, n (x) is unbounded for x > 4 and x < —3. 

13. Additional properties of the Zagier polynomials 

The Zagier polynomials S*(x) have a variety of interesting properties. These 
are recorded here for future studies. 

Coefficients. The Zagier polynomial B*(x) has rational coefficients, some of which 
are integers. Figure [5] shows the number of integer coefficients in B*(x) as a function 
of n. The minimum values seems to occur at the powers 2 J , where the number of 
integer coeffcients is j — 1. 

Signs of coefficients and shifts. The coefficients of B^(x) do not have a fixed 
sign, but there is a tendency towards positivity. Figure 0] shows the excess of 
positive coefficients divided by the total number. On the other hand, the shifted 
polynomial B^(x + |) appears to have only positive coefficients. The coefficients 
of the shifted polynomial appears to be logconcave. This notion is defined in terms 
of the operator L acting on sequences {%■} via C({aj}) — {aj — aj-iaj+i}. A 
sequence is called logconcave if £({oj}) is nonnegative. The sequence is called 
infinitely logconcave if any application of C produces positive sequences. The data 
suggests that the coefficients of B* n {x + |) form an infinitely logconcave sequence. 
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5(1(1 1000 1500 2000 2500 



Figure 2. Integer coefficients 



Figure 3. Linear behavior 



Figure 4. Excess 
of positive coeffi- 
cients 



Figure 5. Roots of B^ooix) 



Roots of £?*. There is a well-established connection between the nature of the 
roots of a polynomials and the logconcavity of its coefficients. P. Branden [3] has 
shown that if a polynomial has only real and negative roots, then its sequence of 
coefficients is infinitely logconcave. This motivated our computations of the roots 



of B*(x). The conclusion is that the polynomial B* (x - 



does not fall in this 



category and Branden's criteria does not apply. Figure [5] shows these roots for 
n = 200. 

A second shift. The polynomial B^(x — |) admits a representation in terms of 
classical special functions. The Gegenbauer polynomial is defined by (see [18], p. 
152, (6.37)): 

(13.1) C^(x)= ( n + 2 n A_1 ) 2 ^ 1 (-n,n + 2A; A+±; 

Theorem 13.1. The shifted Bernoulli polynomial B n (x), defined by B* (x — |) is 

given by 

03*) B„ w =iT„(f) + |:M c a(i). 



Proof. Lemma 19.31 and expanding as a Taylor sum gives 
B n (x) = E 



1 fx l. T 
—Tn — \ — iLn 
n V 2 2 



1 



T, 



fc=i 



iL B 



dx) Tn (2 
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The hypergeometric representation of the Chebyshev polynomial 

(13.3) T n (x) = 2 F 1 (n,-n;±;^) 
and the differentiation rule (Exercise 5.1 in [18], p. 128) 

(13.4) (4-) 2-Fi (a, b; c; z) = ( a ] k ^ k 2 F X (a + k, b + k; c+ k; z) 
give the identity 

(13.5) (J-^ " T n (x) = n2 k ~\k - l)\C™ k (x). 

The odd moments of Lb vanish and the even moments are given by 

B 2k{\) 



. T \ 2k' 

iL B 



<'•>•<» E , 22fe 

according to (Urn]). □ 

(2k) 

The Chebyshev polynomial T n and the Gegenbauer polynomial C n _ 2k have the 
same parity as n. Thus Theorem 113.11 yields a new proof of Theorem lll.il stated 
below in terms of B n . 

Corollary 13.2. The shifted polynomials B n (x) have the same parity as n: 
(13.7) B n {-x) = (-l) n B n (x). 

14. The Euler case 

This section describes a parallel treatment of the Euler polynomial E n (x) defined 
by the generating function 

n— 

Their umbrae is 

2e tx 



(14.2) eval{exp(i<£(x))} = 

The Euler numbers are defined by 

n=0 

and they appear as 

(14.4) E n = 2 n E n (|) . 
Their umbra (£ is 

(14.5) eval(exp(zl£)) = sech(z) 
and the Euler numbers are expressed as 

(14.6) £„ = 2" (<£+!)", 

which is an umbral equivalent of (|14.4|) . 

The next statement is the analogue of Theorem 
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Theorem 14.1. There exists a real valued random variable Le with probability 
density /l e (x) such that, for all admissible functions h, 

(14.7) eval{h(£(x))} = E [h(x - 1/2 + iL E )] 

where the expectation is defined in (|2.12|) . The density of Le is given by 

(14.8) fL E i x ) = sech(7ra;), forxeR. 
In particular, 

(14.9) eval{exp(<£(:c))} = E [it{x - 1/2 + iL E )] 
and 

(14.10) E n {x) =E[(x-± + iL E ) n ] . 

Proof. The proof is similar to the Bernoulli case in Theorem [531 In this case, entry 
3.981.3 of [7J: 



/ sech(ax) cos(xt) dx = — sech f — 
Jo 2a \2a 



(14.11) 
is employed. 

Note 14.2. The analogue of Example 12.51 is 

(14.12) eval{log<£(a;)} = log 2 + 2 logT (^) -21o g r(§) 

and differentiation produces 



□ 



for x > h 



(14.13) 

with 

(14.14) 



eval{£- fe (x)} = LJ^2^ k ~ 1 \x). 



(3(x) 



(k-iy 



for x > 



the beta function on page 906 of [7J. The proofs of all these results are similar to 
those presented for the Bernoulli case. 



It is natural to consider now the modified Euler numbers 

n 



(14.15) 



r=0 



n + r 
2r 



-EL 



n > 0. 



n 



Symbolic experimentation suggested the next statement. The proof of the next 
statement follows the same ideas as in the Bernoulli case. 

Theorem 14.3. The odd subsequence of the modified Euler numbers {-Eln+i} * s a 
periodic sequence of period 3 ; with values {1, —2, 1}. 

Define the modified Euler polynomials by 



(14.16) 



r=0 



n + r 



2r / n + r 



E r (x). 



Then the even-order numbers -Ej^O) have period 12 with values 



n mod 12 2 


4 


6 


8 


10 




i 





-2 


3 


-2 
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The proof of this result follows the same steps as in the Bernoulli case. 

The final statement in this section is the analogue of Theorem 111.11 
Theorem 14.4. The modified Euler polynomials satisfy 
(14.17) E* n (-x-3) = (-irE* n (x). 

15. The duplication formula for Zagier polynomials 
The identity 

(15.1) B k (mx) ^ m k - 1 B k [ x + — 

was given by J. L. Raabe in 1851. The special case m = 2 gives the duplication 
formula for Bernoulli polynomials 

(15.2) 2B k (2x) = 2 k B k (x) + 2 k B k (x + ±) . 

Summing over k yields 
(15.3) 

A fn + k\ B k (2x) _ -A /n + fc\ 2 fc B fc (x) A /n + fc\ 2 fc 5 fc (a; + \) 
2-^\2k)n + k~2-^\2k)n + k + 2-^\2k) n + k 

An umbral interpretation of this identity leads to a duplication formula for the 
Zagier polynomials. This result is expressed in terms of the umbral composition 
defined next. 

Definition 15.1. Given two sequences of polynomials P = {P n (x)} and Q = 
{Qn(x)}, their umbral composition is defined as 

n 

(15.4) {P°Q)n{x) = J2Pk,nQk(x), 

k=0 

where Pk,n is the coefficient of x k in P n {x). 

The use of umbral composition is clarified in the next lemma. 
Lemma 15.2. Let P and Q be polynomials and assume 

(15.5) P n (x) = eval{(x + and Q n {x) = eval{ (a; + 0) n } . 
Then 

(15.6) (PoQ) n (x) =eval{(a; + f P + Q)™}. 
Proof. Denoting the relevant umbrae by a subindex, then 

evabp^Uz + fp + Q)™} = eval Q {P„(x + 0)} 

71 

= ^Pk,nQk(x) 
fc=0 

= (PoQ) n (x), 

as claimed. □ 
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Consider now the Bernoulli and Euler umbrae 

t 2 

(15.7) eval{exp(i<8)} = — and eval {exp(i£)} = — - 

e 1 — 1 e J + 1 

given in (|2.7[) and (|14.2I) . respectively. The identity 

(15.8) eval {exp(£*8)} x eval {exp(i£)} = eval {exp(2i<8)} 
is written (at the umbrae level) as 

(15.9) <B + £ = 2<B. 

The first summand on the right of (|15.3[) contains the term 
2 k B k (x) = eval{2 fe (:r + <B) fe } 
= eval{(2a; + 2<B) fc } 
= eval{(2x + < B + £) fe } 
= eval{(«Bo (g) fc (2x)}. 
Lemma 115.21 has been used in the last step. Similarly 

2 k B k {x + \)= cval{(<B o €) k (2x + 1)} . 

Thus, (fTOl) reads 

(15.10) 2K( 2x ) = ( B * ° E)„(2x) + (B* o E) n (2x + 1) 
that can also be expressed in the form 

(15.11) 2B; i {2x) = (B* o E(x)) n {x) + (B* o E (x + |)) n (x + ±) , 
that is an analogue of (|15.2[) for the Zagier polynomials. 
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